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We present the application of the variational-wavelet analysis to the quasiclassical 
calculations of the solutions of Wigner/von Neumann/Moyal and related equa- 
tions corresponding to the nonlinear (polynomial) dynamical problems. (Naive) 
deformation quantization, the multircsolution representations and the variational 
approach are the key points. We construct the solutions via the multiscale ex- 
pansions in the generalized coherent states or high-localized nonlinear eigenmodes 
in the base of the compactly supported wavelets and the wavelet packets. We 
demonstrate the appearance of (stable) localized patterns (waveletons) and con- 
sider entanglement and decoherence as possible applications. 



1. Wigner-like Equations 

In this paper we consider the applications of a numerical-analytical tech- 
nique based on local nonlinear harmonic analysis (wavelet analysis, general- 
ized coherent states analysis) to the quasiclassical calculations in nonlinear 
(polynomial) dynamical problems in the Wigncr-Moyal approach. The cor- 
responding class of Hamiltonians has the form 

H(p,q) = ^ + U(p,q), (1) 

where U (p, q) is an arbitrary polynomial function on p, q, and plays the 
key role in many areas of physics [1], [2]. The particular cases, related to 
some physics models, are considered in [3]-[12]. Our goals are some attempt 
of classification and the explicit numerical-analytical constructions of the 
existing quantum states in the class of models under consideration. There 
is a hope on the understanding of relation between the structure of initial 
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Hamiltonians and the possible types of quantum states and the qualita- 
tive type of their behaviour. Inside the full spectrum there are at least 
three possibilities which arc the most important from our point of view: 
localized states, chaotic-like or/and entangled patterns, localized (stable) 
patterns (definitions can be found below). All such states are interesting in 
the different areas of physics (e.g., [1], [2]) discussed below. Our starting 
point is the general point of view of a deformation quantization approach 
at least on the naive Moyal/Weyl/Wigner level [1], [2]. The main point of 
such approach is based on ideas from [1] , which allow to consider the alge- 
bras of quantum observables as the deformations of commutative algebras 
of classical observables (functions). So, if we have as a model for classical 
dynamics the classical counterpart of Hamiltonian (1) and the Poisson man- 
ifold M (or symplectic manifold or Lie coalgebra, etc) as the corresponding 
phase space, then for quantum calculations we need first of all to find an 
associative (but non-commutative) star product * on the space of formal 
power series in h with coefficients in the space of smooth functions on M 
such that 

f*g = fg + Hf, g} + J2 s), (2) 

n>2 

where {/, g} is the Poisson brackets, B n are bidiffcrcntial operators. Kont- 
sevich gave the solution to this deformation problem in terms of the formal 
power series via the sum over graphs and proved that for every Poisson 
manifold M there is a canonically defined gauge equivalence class of star- 
products on M. Also there are the nonperturbative corrections to power 
series representation for * [1]. In the naive calculations we may use the 
simple formal rule: 

* EECX p(y( ,9 <J ,9 P~ dpdqi) ( 3 ) 

In this paper we consider the calculations of the Wigner functions W(p, q, t) 
(WF) corresponding to the classical polynomial Hamiltonian H(p,q,t) as 
the solution of the Wigner- von Neumann equation [2]: 

ih— W = H * W - W * H (4) 

and related Wigner-like equations. According to the Weyl transform, a 
quantum state (wave function or density operator p) corresponds to the 
Wigner function, which is the analogue in some sense of classical phase- 
space distribution [2]. We consider the following form of differential equa- 
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tions for time-dependent WF, W — W(p, q, t): 



Wt = — sin 



■HW 



'i<d H d w -d H d w ) 

2 \Vq "p w p l 

which is a result of the Weyl transform of von Neumann equation 

dp 



ih 



dt 



[H, P ] 



(5) 



(6) 



In our case (1) we have the following decomposition [2] (U = U(q) in the 
following only for simplicity, but the case U = U(p, q) can be considered 
analogously): 



where 



u 



(21 + 1)! dq 2 ^ 1 d v 2l+x 



(7) 



(8) 



-P-^-W(p,q;t) 
m oq 



(9) 



Let {\E),E} be the full set of discrete/continuous eigenfunctions (eigenval- 
ues) 



H\E) = E\E) 

then we have the following representation for the Moyal function: 



(10) 



W E ".E'(q,p) 



1 

2^h 



1 



dZe- i '*' h (q+-Z\E"){E'\q--Z) (11) 



which is reduced to the standard WF in the case E' = E": We,e(p,q) = 
W(p,q). As a result, the time independent Moyal function generates the 
time evolution of the WF. The corresponding integral representation con- 
tains the initial value of the density operator p(0) as a factor [2]. The Moyal 
function satisfies the following system of (pseudo) differential equations 



2m 8m dq 2 



£ 

i=i 



(-iy(h/2) 21 d 2l u & 



[2l)\ dq 21 dp 21 



WE- 



S' + E" 



We" 



(12) 
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_dUd_ _Y^ (-i) l (h/2) 21 d 2l+1 U d 2l+1 

in fin An r)n — ^ 



mdq dqdp (21 + 1)! dq 2l + x dp 2l+1 



W E » 



= 1 -{E" -E')W e »,e> (13) 

really nonlocal/pseudodiffcrential for arbitrary Hamiltonians. But in case 
of polynomial Hamiltonians (1) we have only a finite number of terms in 
the corresponding series. Also, in the stationary case after Weyl-Wigner 
mapping we have the following equation on WF in c- numbers [2]: 

( p 2 h p 8 h 2 8 2 \ . , 

fc + 2^-^r M + (14) 

il d )w(p,q)=eW( P ,q) 



U(q 



2i dp> 

After expanding the potential U into the Taylor series we have two real 
partial differential equations which correspond to the mentioned before par- 
ticular case of the Moyal equations (12), (13). 

Our approach, presented below, in some sense is motivated by the analy- 
sis of the following standard simple model considered in [2] . Let us consider 
model of interaction of nonresonant atom with quantized electromagnetic 
field: 

v 2 

H=-^ + U(x), U{x) = U {z,t)g(x)a+a (15) 

where potential U depends on creation/annihilation operators and some 
polynomial on x operator function (or approximation) g(x). It is possible 
to solve Schroedinger equation 

d|* > 

ih^-— = jt*> (16) 

by the simple ansatz 

|*(i) >= ^2 Wn / dx \^n{x, t)\x > ®|n > (17) 

-oo 

which leads to the hierarchy of analogous equations with potentials created 
by n-particle Fock subspaces 



ih 



-2 

\^ + U (t)g(x)nU n (x,t) (18) 
1 2m J 



dt 

where vE , n (x, t) is the probability amplitude of finding the atom at the time 
t at the position x and the field in the n Fock state. Instead of this, we 
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may apply the Wigner approach starting with proper full density matrix 

p=|*(*) ><*(*)! = (19) 
Wn'Wn" j dx ' J dx"^ n/ (x',t)^* n „(x",t)\x' >< x"\ ® \n' >< n"\ 

n' ,n" 

Standard reduction gives pure atomic density matrix 

Pa= < n\p\n > = (20) 

Jn=0 

J2\ W n\ 2 j Ax' J dx"y n (x\t)y* n (x",t)\ X ' ><X>'\ 

Then we have incoherent superposition 

oo 

W(x,p,t) = Y / \w n \ 2 W n {x,p,t) ( 21 ) 

n=0 

of the atomic Wigner functions 

W n (x,p,t) = ^ J dtexp(-^)** n (x-y,t)* n (x+^,t) (22) 

corresponding to the atom motion in the potential U n {x) (which is not more 
than polynomial in x) generated by n-level Fock state. They are solutions 
of proper Wigner equations: 

dWn V dWn ^ (-l)'(fi/2) 2 ' d 2i+1 U n (x) d 2l!+1 W n 
8t m dx f^ Q (2£+l)\ dx 2 ^ 1 d P 2t+1 [ ' 

In the following section we'll generalize this construction and we are 
interested in description of entanglement in the Wigner formalism. 

The next example describes the decoherence process. Let we have col- 
lective and environment subsystem with their own Hilbert spaces 

n = n c ®n e (24) 

Relevant dynamics are described by three parts including interaction 

H = H c ®I e +I c ®H e + H int (25) 
For analysis, we can choose Lindblad master equation [2] 

p=^ h [H,p}-J27n(L+L n p + pL+L n -2L nP L+) (26) 



February 1, 2008 6:42 WSPC/Trim Size: 9in x 6in for Proceedings 



fz 



7 



which preserves the positivity of density matrix and it is Markovian but 
it is not general form of exact master equation. Other choice is Wigner 
transform of master equation [2] and it is more preferable for us 



In the next section we consider the variation-wavelet approach for the 
solution of all these Wigner-like equations (6)-(9), (12), (14), (23), (26), 
(27) for the case of an arbitrary polynomial U(q,p), which corresponds to a 
finite number of terms in the series from (8), (12), (13), (14), (23), (27) or 
to proper finite order of h. Our approach is based on the extension of our 
variational- wavelet approach [3]- [12]. Wavelet analysis is some set of math- 
ematical methods, which gives the possibility to work with well-localized 
bases in functional spaces and gives maximum sparse forms for the general 
type of operators (differential, integral, pseudodifferential) in such bases. 
These bases are the natural generalization of standard coherent, squeezed, 
thermal squeezed states [2] , which correspond to quadratical systems (pure 
linear dynamics) with Gaussian Wigner functions. Because the affine group 
of translations and dilations (or more general group, which acts on the space 
of solutions) is inside the approach (in wavelet case), this method resembles 
the action of a microscope. We have a contribution to the final result from 
each scale of resolution from the whole underlying infinite scale of spaces. 
Our main goals are an attempt of classification and construction of possible 
nontrivial states in the system under consideration. We are interested in the 
following states: localized, entangled patterns, localized (stable) patterns. 
We start from the corresponding definitions (at this stage these definitions 
have only qualitative character). 

1. By localized state (localized mode) we mean the corresponding (par- 
ticular) solution of the system under consideration which is localized in 
maximally small region of the phase space. 

2. By chaotic/entangled pattern we mean some solution (or asymptotics 
of solution) of the system under consideration which has equidistribution 
of energy spectrum in a full domain of definition. 

3. By localized pattern (waveleton) we mean (asymptotically) stable 
solution localized in relatively small region of the whole phase space (or a 
domain of definition). In this case all energy is distributed during some 
time (sufficiently large) between few localized modes (from point 1) only. 

Numerical calculations explicitly demonstrate the quantum interference 



W = {H,W} PB + 



(27) 




dl n+1 U(x)d 2 p n+1 W{x, P ) + 2jd p pW + Dd 2 p W 
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of generalized coherent states, pattern formation from localized eigenmodes 
and the appearance of (stable) localized patterns (waveletons). 

2. Variational Multiscale Representation 

We obtain our multiscale/multiresolution representations for solutions of 
Wigner-like equations via a variational-wavelet approach. We represent 
the solutions as decomposition into modes related to the hidden underlying 
set of scales [13]: 

oo 

W(t,q,p) = @5 i W(t,q,p) (28) 

where value i c corresponds to the coarsest level of resolution c or to the 
internal scale with the number c in the full multiresolution decomposition 
of underlying functional space (L 2 , e.g.) corresponding to problem under 
consideration: 

V c C V c+1 C V c+2 C . . . (29) 

andp = (pi,p 2) -),q= (<?i,<72, ■ ••)>£» = (pi,<?i, -;Pi,qi) arc coordinates in 
phase space. In the following we may consider as fixed as variable numbers 
of particles. The second case corresponds to quantum statistical ensemble 
(via "wignerization" procedure) and will be considered in details elsewhere 
[12]. 

We introduce the Fock-like space structure 

^ = E® ff " ( 3 °) 

for the set of n-particle wave functions (states): 

W = {W , W!{ Xl ;t), W 2 (x 1 ,x 2 ;t), W N {x u . . . ,x N ;t), . . .}, (31) 

where W p (xi, . . . , x p ; t) G HP, H° = C, H p = L 2 (R 6 p) (or any different 
proper functional space), W G H°° = H° H 1 © • • • H p © . . . with 
the natural Fock space like norm (guaranteeing the positivity of the full 
measure) : 

(W,W) = W 2 + Y, f W?(x u ...,x i ;t)f[M. (32) 

First of all we consider W = W(t) as a function of time only, W £ L 2 {R), 
via multiresolution decomposition which naturally and efficiently introduces 
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the infinite sequence of the underlying hidden scales [13]. We have the con- 
tribution to the final result from each scale of resolution from the whole infi- 
nite scale of spaces. We consider a multiresolution decomposition of L 2 (R) 
(of course, we may consider any different and proper for some particular 
case functional space) which is a sequence of increasing closed subspaces 
Vj E L 2 (R) (subspaces for modes with fixed dilation value): 

...V-2 C V-i C V C Vi c V 2 C ... (33) 

The closed subspace Vj(j e Z) corresponds to the level j of resolution, or to 
the scale j and satisfies the following properties: let Wj be the orthonormal 
complement of Vj with respect to Vj + \\ Vj + \ = Vj Wj. Then we have 
the following decomposition: 



{W(t)}= Wj or {W(t)} = V ($W j , (34) 

— oo<j<oo j— 

in case when Vq is the coarsest scale of resolution. The subgroup of trans- 
lations generates a basis for the fixed scale number: spa,n kGZ {2^ 2 ^(2H — 
k)} = Wj. The whole basis is generated by action of the full affine group: 

span keZjeZ {V^(2H - k)} = span fciJ . eZ {* J -, fc } - {W(t)}. (35) 

Let the sequence {Vj}, Vj C L 2 (R) correspond to multiresolution analysis 
on the time axis, {V Xi } correspond to multiresolution analysis for coor- 
dinate Xi, then Vj 1 ^ = Vf 1 (g> • • • <g> Vf n <g> Vj corresponds to the mul- 
tiresolution analysis for the n-particle function W n (x\, . . . , x n ; t). E.g., 
for n = 2: Vq = {/ : f(x 1 ,x 2 ) = Efei,fe 2 ^i.fe^^i - fc i^2 - 
k 2 ), a klM e £ 2 (Z 2 )}, where (f> 2 (x 1 ,x 2 ) = (j) 1 {xi)<t> 2 {x 2 ) = (j) 1 ® 4> 2 {xi,x 2 ), 
and <j) l {xi) = 4>{xi) form a multiresolution basis corresponding to {Vp}. 
If {</> 1 (xi — £)}, i 6 Z form an orthonormal set, then <j) 2 {x\ — ki,x 2 — k 2 ) 
form an orthonormal basis for V 2 . So, the action of the affine group gener- 
ates multiresolution representation of L 2 (R 2 ). After introducing the detail 
spaces W 2 , we have, e.g. V 2 = V 2 © Wq . Then the 3-component basis for 
Wq is generated by the translations of three functions 

*? = <j>\ Xl ) ® <f 2 (x 2 ), tf| = 2 (a , 2)j ^2 = ^ &( X2 y 

Also, we may use the rectangle lattice of scales and one-dimensional wavelet 
decomposition: 

f{x u x 2 ) = {f^i,i®^i,k)^i,e®^j,k{xi,x 2 ), 

i,£;j,k 
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where the basis functions "J/^ <8> k depend on two scales 2~ l and 2~ J . 
After constructing the multidimensional basis we may apply one of the 
variational procedures from [3]-[12]. We obtain our multiscale/multireso- 
lution representations (formulae (40) below) via the variational wavelet 
approach for the following formal representation of the systems from the 
Section 1 (or its approximations). 

Let L be an arbitrary (non)linear differential/integral operator with 
matrix dimension d (finite or infinite), which acts on some set of func- 
tions from L 2 (Vl® n ): * = x\, x 2 , ■ ■ ■ ) = Xi, x 2 , ...),... , 

^ d (t,Xi,X2, •••))> Xi e £1 (1 R 6 , n is the number of particles: 

L* = L(Q,t,Xi)*(t,Xi)=0, (36) 

where 

Q = Qd ,di,d2,-{^ x ux 2 , • • • ,d/dt, d/dxx, d/dx 2 , . . . , J Hk) = 

fd\v>( d d \i2 f 

Let us consider now the N mode approximation for the solution as the 
following ansatz: 

N 

V N (t,x 1 ,x 2 ,...)= <k i 1 i2...A io ®B il ®Ci 2 ...(t,x 1 ,X2,...). 

(38) 

We shall determine the expansion coefficients from the following conditions 
(different related variational approaches are considered in [3]- [12]: 



Cfci,fc 2 ,...= J (L^)A ka (t)B kl (x 1 )C k2 (x 2 )dtdx 1 dx 2 --- = 0. (39) 

Thus, we have exactly dN n algebraical equations for dN n unknowns 
fflio,»i,...- This variational approach reduces the initial problem to the prob- 
lem of solution of functional equations at the first stage and some algebraical 
problems at the second. We consider the multiresolution expansion as the 
second main part of our construction. So, the solution is parametrized 
by the solutions of two sets of reduced algebraical problems, one is linear 
or nonlinear (depending on the structure of the operator L) and the rest 
are linear problems related to the computation of the coefficients of the 
algebraic equations (39). These coefficients can be found by some wavelet 
methods by using the compactly supported wavelet basis functions for the 
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expansions (38). As a result the solution of the equations from Section 1 
has the following multiscale or multiresolution decomposition via nonlinear 
high-localized eigenmodes 

W(t,X!,X 2 , ■ ■ ■) = ^2 ClijU 1 ®V 3 (t,Xi,X2,..-), 

V 3 (t) = V^ low (t) + Y,V l i ^it), wi~2', (40) 
;>jv 

U\x s ) = U l M slow (x s ) + ]T lP m (k' m x a ), k s m ~ 2™, 

m>M 

which corresponds to the full multiresolution expansion in all underlying 
time/space scales. The formulae (40) give the expansion into a slow part 
and fast oscillating parts for arbitrary N, M. So, we may move from the 
coarse scales of resolution to the finest ones for obtaining more detailed 
information about the dynamical process. In this way one obtains contri- 
butions to the full solution from each scale of resolution or each time/space 
scale or from each nonlinear eigenmode. It should be noted that such repre- 
sentations give the best possible localization properties in the correspond- 
ing (phase)space/time coordinates. Formulae (40) do not use perturbation 
techniques or linearization procedures. Numerical calculations are based 
on compactly supported wavelets and related wavelet families [13] and on 
evaluation of the accuracy on the level N of the corresponding cut-off of 
the full system regarding norm (32): 

\\W N+1 -W N \\<e. (41) 

So, by using wavelet bases with their best (phase) space/time local- 
ization properties we can describe localized (coherent) structures in quan- 
tum systems with complicated behaviour. The modeling demonstrates the 
appearance of different (stable) pattern formation from high-localized co- 
herent structures or chaotic behaviour. Our (nonlinear) eigenmodes are 
more realistic for the modelling of nonlinear classical/quantum dynamical 
process than the corresponding linear gaussian-like coherent states. Here 
we mention only the best convergence properties of the expansions based 
on wavelet packets, which realize the minimal Shannon entropy property 
and the exponential control of convergence of expansions like (40) based 
on the norm (32). Fig. 1 shows the high- localized eigenmode contribu- 
tion to the WF, while Fig. 2, 3 give the representations for the full solu- 
tions, constructed from the first 6 eigenmodes (6 levels in formula (40)), 
and demonstrate the stable localized pattern formation (wavcleton) and 
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complex chaotic-like behaviour. Fig. 3 corresponds to (possible) result of 
superselection (einselection) [2] after decoherence process started from Fig. 
2 or Fig. 4. Fig. 5 and Fig. 6 demonstrate time steps during appearance 
of entangled states. It should be noted that we can control the type of 
behaviour on the level of the reduced algebraical system (39) [12]. 
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Figure 1. Localized mode contribution to Wigner function. 
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Figure 2. Chaotic-like Wigner function. 
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Figure 3. Localized pattern-like (waveleton) Wigner function. 
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Figure 4. Entangled-like Wigner function. 
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Figure 5. Section of the coarse level approximation for Wigner function. 
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Figure 6. 



Section of the finest level approximation for Wigner function. 



